Mon. Not. R. Astron. Soc. 000, 1-10 (200X) Printed 28 October 2008 (MN HT E X style file vl.4) 



00 
O 

o 

(N 



The radio remnant of SN1993J: an instrumental 
explanation for the evolving complex structure 

Ian Hey wood 1 , Katherine M. Blundell 1 , Hans-Rainer Klockner 1 , Anthony J. Beasley z 

1 University of Oxford Astrophysics, Keble Road, Oxford, 0X1 3RH, UK 

2 NEON, Inc., 3223 Arapahoe Ave., Suite 210, Boulder, CO 80303, USA 

Accepted 2008 month 10. Received 200X month XX. Preprint version for astro-ph 

ABSTRACT 

We present simulated images of Supernova 1993J at 8.4 GHz using Very Long 
Baseline Interferometry (VLBI) techniques. A spherically symmetric source model is 
convolved with realistic uu-plane distributions, together with standard imaging proce- 
dures, to assess the extent of instrumental effects on the recovered brightness distribu- 
tion. In order to facilitate direct comparisons between the simulations and published 
VLBI images of SN1993J, the observed uv-coverage is determined from actual VLBI 
observations made in the years following its discovery. 

The underlying source model only exhibits radial variation in its density profile, 
with no azimuthal dependence and, even though this model is morphologically sim- 
ple, the simulated VLBI observations qualitatively reproduce many of the azimuthal 
features of the reported VLBI observations, such as appearance and evolution of com- 
plex azimuthal structure and apparent rotation of the shell. We demonstrate that such 
features are inexorably coupled to the uu-plane sampling. 

The brightness contrast between the peaks and the surrounding shell material are 
not as prominent in the simulations (which of course assume no antenna- or baseline- 
based amplitude or phase errors, meaning no self-calibration procedures will have 
incorporated any such features in models). It is conclusive that incomplete uu-planc 
sampling has a drastic effect on the final images for observations of this nature. Dif- 
ference imaging reveals residual emission up to the 8a level. Extreme care should be 
taken when using interferometric observations to directly infer the structure of objects 
such as supernovae. 

Key words: stars: supernovae - supernovae: individual (SN1993 J) - radio continuum: 
stars - techniques: interferometry 



1 INTRODUCTION 

Supernova 

1993J (RA = 09h 55m 24.7747s, 8 = +69° 01' 13.7031") 
was discovered on 28 March 1993 in the nearby spiral galaxy 
M81 (Ripero et al., 1993) and became the brightest radio su- 
pernova ever observed. Its brightness and high Declination 
have allowed this object to be studied in great detail at radio 
wavelengths. 

A shell-like radio structure associated with SN1993J 
was resolved using VLBI techniques as early as 239 days 
after the supernova explosion (Marcaide et al., 1995). An 
extensive VLBI monitoring program then followed which 
tracked the expansion of the shell over seven years at three 
frequencies (Bietenholz, Bartel & Rupen, 2001; Bartel et 
al., 2002; Bietenholz, Bartel & Rupen, 2003). These obser- 
vations yielded a wealth of information about the supernova, 



including accurate determination of the explosion centre, 
measurement of the deceleration of the ejecta and radio light 
curves and spectra. Images of the resolved shell allowed di- 
rect monitoring of its development, from its initial barely 
resolved state, through expansion into partial shell struc- 
tures, to the complex structure apparently present in the 
final epochs 3000+ days later. 

Behaviour such as apparent rotation of features in the 
shell in the plane of the sky are seen between epochs of VLBI 
observations, as are phenomena such as coherent structures 
dividing into a series of peaks which apparently change in 
azimuthal position as the remnant expands. Bietenholz, Bar- 
tel and Rupen (2003) discuss possible physical mechanisms 
which could explain the structure in the shell, including 
Rayleigh- Taylor cooling instabilities, evolving magnetic field 
structure and interactions with inhomogeneous circumstel- 
lar material. 
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This paper presents a new investigation of this intrigu- 
ing evolution of the brightness distribution in the expand- 
ing shell. We explore whether the development of apparently 
complex structure may be due to incomplete sampling of the 
uv plane, which is problematic especially for sparse interfer- 
ometer arrays such as VLBI networks. Similar effects had 
also been a cause for concern when interpreting 6-cm MER- 
LIN images of nova Cassiopeiae 1995 (Heywood et al. 2005; 
Heywood & O'Brien, 2007) where the apparently complex 
development of the radio emission could be accounted for by 
a spherically symmetric shell model and sparse sampling of 
the mi-plane. With this in mind wc simulated 8.4-GHz VLBI 
observations of the radio shell of SN1993J at nine epochs, 
for comparison with the genuine VLBI images. 

Further discussion regarding the non-uniqueness of de- 
convolved radio images (particularly due to the CLEAN de- 
convolution method) may be found in e.g. Cornwell, Braun 
and Briggs (1999) and Briggs (1995). The latter contains 
a particularly relevant case study of the radio shell of 
SN1987A whereby different deconvolution methods yield 
subtle variations in apparent structure. 



2 SIMULATIONS 

Creation of the simulated VLBI images requires two com- 
ponents. First, a model of the assumed source brightness 
distribution is required. This model should represent how 
an ideal instrument would image the source at a given fre- 
quency. The second component of the simulation is a uv- 
plane sampling function which accurately represents that of 
the real instrument. 

2.1 Model brightness distributions 

A synchrotron shell model is used as the basis for our simu- 
lated radio images. The expansion of the shell in this model 
is not self-similar because the radial density profile is time 
dependent as the physical shell boundary expands, and re- 
verse shocks propagate back through the ejecta. The model 
is fully spherically symmetric however, and the density 
variations are purely radial with no azimuthal component. 
A detailed description of this model is provided by Mio- 
duszewski et al. (2001). Simulated radio images were used 
with ages that closely matched the epochs of the VLBI ob- 
servations, and were evenly spaced throughout the evolution 
of SN1993J. Mioduszewski et al. (2001) also parametrized 
and fitted their models to the radio observations to give 
consistent angular sizes and flux density values, and kindly 
allowed us to use their models in the investigations presented 
here. Figure 1 shows three epochs of the simulated radio im- 
ages and their corresponding flux density profiles. 

2.2 Recreating the uv coverage 

In addition to the model images being chosen to coincide 
as best as possible with the genuine VLBI observations the 
model array for each epoch is defined to be the same as 
those used during the SN1993J VLBI campaign. Only sta- 
tions which were used during an actual observation are in- 
cluded for the corresponding simulation. 

The AIPS task UVCON is used to generate the model 



tro-data sets. By examining VLBA archive files, the start and 
end time for any observation are used to determine the total 
duration. Using the relevant model image and an antenna 
file describing the stations present for a specific observation, 
UVCON is executed, resulting in a track length equivalent 
to this total duration. In the case of the real observation, 
this track would encompass all sources and frequencies. 

In order to reduce the u-u-data set such that it only 
includes visibilities related to SN1993J at 8.4 GHz three 
flag tables are applied, documenting the frequency switching 
information, the SN1993J scan times, and the automated 
flagging that occurred at the receivers. 

These flag tables are constructed from archival data 
files, and converted to a format that can be directly read 
by the UVFLG task in AIPS. Bad data that are manu- 
ally removed from the true VLBI observations in the pre- 
calibration stage cannot be accounted for as the simulations 
are not prone to problems such as radio frequency inter- 
ference. This means that in each case the simulations are 
over-estimating the uw-plane sampling, and that our analy- 
sis errs on the side of being conservative. 

Multi-frequency interferometric observations often em- 
ploy a technique to improve uw-plane coverage, whereby the 
observing frequency is switched periodically. Consider an 
18-hour track at three frequencies. Rather than observing 
at each frequency in three successive six-hour blocks the 
observing frequency is switched, for example every 20 min- 
utes. The visibilities for a given frequency are therefore more 
spread out on the ttv-plane, facilitating the sampling of more 
diverse Fourier components and reducing linear biases. De- 
termining this frequency switching information allows us to 
flag out scans which do not occur at 8.4 GHz. 

Once the simulated uw-tracks are consistent only with 
the scans at 8.4 GHz we flag out times where the calibration 
sources are being observed. From the archive data we can 
obtain the times when SN1993J was being observed and use 
this information to retain only these time ranges in the sim- 
ulation. AlPS-readable flag tables describing data that were 
flagged at the time of the observation (e.g. due to slewing or 
pointing errors, source changes etc.) are also available, and 
this is the final flag table applied to the simulated data set. 

Note also that a 15° elevation cut off for each antenna is 
also imposed on the simulation. Data are routinely flagged 
below this elevation because system temperatures rise as the 
antenna points closer to the horizon, and local geographical 
features can often obscure the line of sight. For a source 
at the Declination of SN1993J this condition only slightly 
affects the longest baselines in the VLBA network. 

The UVCON task also adds random Gaussian noise to 
the simulated data set, the level of which is governed by the 
parameters of each antenna in the array, such as the antenna 
efficiency and system temperature. 

2.3 Cleaning and contouring 

The simulated uw-data sets were imaged using the CLEAN 
algorithm (Clark, 1980) using the AIPS task IMAGR. For 
each simulation the circular restoring beam used is consis- 
tent with that of the nearest epoch of the VLBI observations, 
as reported by Bietenholz, Bartel and Rupen (2003). Briggs 
(1995) weighting was utilised with the ROBUST parameter 
in IMAGR set to zero. For the nine simulated maps pre- 
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Figure 1. Model radio brightness distributions (upper row) and their corresponding brightness profiles (lower row) for three epochs of 
the parametrized synchrotron model for SN1993J (Mioduszewski et al., 2001). These models are not self-similar, having radial evolution, 
but they are entirely azimuthally symmetric. 



sented in Figure 2 the data were reweighted by the square 
roots of the data weights using the UVWTFN parameter. 
This is a common technique for VLBI observations and is 
consistent with the imaging methods of Bietenholz, Bartel 
and Rupen (2003). Deconvolution was carried out interac- 
tively and the clean cycle terminated when it was judged 
that all substantial flux has been removed. 



3 SIMULATED OBSERVATIONS 

In this section we first present the nine simulated radio im- 
ages which best recreate the genuine VLBI observations of 
SN1993J. We discuss aspects of the real VLBI observations 
which differ from our simulations. Also presented in this 
section are residual images created by differencing the sim- 
ulated radio image from a convolved input model. 



Since no antenna-based or baseline-based amplitude or 
phase errors are included in our simulations, this means no 
self-calibration procedures will be responsible for spuriously 
enhancing features in our simulations. Our simulations will 
only underestimate the extent to which small compact fea- 
tures could become reinforced in self-calibration and clean- 
ing cycles. 

Contouring of the simulated images is also carried out in 
a way that is consistent with the images presented by Bieten- 
holz, Bartel and Rupen (2003). The greyscale in each image 
runs from 8% to 100% of the peak brightness for epochs 
below 1253 days and from 16% to 100% for epochs above. 
Contour levels are 1%, 2%, 4%, 8%, 16%, 32%, 45%, 64% and 
90% of the peak brightness in each case, starting with the 
first contour level that is brighter than three times the root 
mean square background brightness levels (3a). Presenting 
the simulated images in a consistent way better facilitates 
comparison with the real observations. 



3.1 Simulated VLBI Images 

Figure 2 shows the resulting images from nine simulated 
VLBI observations of SN1993J. A summary of the relevant 
image parameters can be found in Table 1. The total shell 
flux density values for the simulated VLBI maps are mea- 
sured by using the AIPS task TVSTAT to define the region 
of emission. 

3.2 Evolution of the simulated radio images 

The evolution of the simulated VLBI images as shown in Fig- 
ure 2 has much in common with that of the genuine VLBI 
observations of Bietenholz, Bartel and Rupen (2003), with 
the major difference occuring in the early stages. Specifically, 
in both the case of the real observations and the simulations 
the emission is intially barely resolved, and as the shell ex- 
pands the real observations begin to take on a partial shell 
structure, which is only hinted at in the simulated image 
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Table 1. Summary of the image parameters for the input model and the simulations. 



Observation 


Input model 




Simulated VLBI images 




Restoring 


time (days) 


Time (days) 


Flux density (mjy) 


Flux density (: 


mjy) Peak (mjy/beam) 


3<r (mjy/beam) 


beam (mas) 


352 


321 


86.28 


85.56 


12.24 


0.087 


0.57 


390 


390 


71.86 


70.74 


9.21 


0.143 


0.60 


686 


707 


38.79 


38.34 


3.43 


0.038 


0.80 


996 


867 


31.06 


30.82 


3.04 


0.213 


0.93 


1107 


1062 


25.57 


24.86 


1.99 


0.148 


0.97 


1532 


1539 


19.26 


17.01 


1.56 


0.276 


1.10 


1893 


1939 


18.32 


17.91 


0.99 


0.066 


1.12 


2525 


2351 


16.84 


16.38 


0.68 


0.037 


1.12 


2787 


2846 


15.15 


14.81 


0.61 


0.067 


1.12 



corresponding to day 996. While in the VLBI observations 
this partial shell structure expands and appears to rotate, 
the simulated images deviate very little from circularity un- 
til day 1532 when a fragmented shell develops. 

The fragmented shell continues to expand, giving the 
impression of ridge structures dividing into peaks and un- 
dergoing apparent rotation. The simulation with the largest 
angular extent is remarkably similar to the VLBI observa- 
tion at day 2787. 

Although the ring-like structure upon which the peaks 
are situated seems far more coherent in the simulations they 
are still qualitatively very similar to the VLBI images. Given 
that the highest contour represents 90% of the peak bright- 
ness and the next lowest one is at the 64% level these bright- 
ness enhancements represent a significant deviation from cir- 
cularity, which given the nature of our input model can only 
be a purely instrumental effect. If the structure apparent 
in the VLBI images of SN1993J represents genuine dumpi- 
ness of the ejecta then at best a degeneracy exists between 
structures of that nature and ring-like objects. 

3.3 How the model deviates from the reality 

Our simulated images can qualitatively reproduce the ef- 
fects apparent in the VLBI observations of SN1993J, how- 
ever some differences are present and this section suggests 
some possible reasons why such differences may arise. Devi- 
ations from a perfect reconstruction of the SN1993J VLBI 
observations will occur either in the uw-plane model or the 
radio brightness distribution model. 

An obvious reason why the uw-distributions will devi- 
ate is that our simulations are not prone to radio frequency 
interference. Any contamination of the real data by such in- 
teference would be flagged out, therefore in this respect our 
simulations are over-estimating the sampling of the uv plane. 
In addition to this, as mentioned elsewhere, the simulations 
are not prone to calibration errors which may result in spu- 
rious features that are reinforced due to self-calibration. 

With respect to the model brightness distribrutions em- 
ployed, it may well be the case that the radio remnant 
of SN1993J is not spherically symmetric with solely ra- 
dial structure variations. Indeed, imaging nearby supernova 
remnants at many wavelengths suggests that complex mor- 
phologies are common. However, should it be the case that 
SN1993J deviates from sphericity as much as the VLBI im- 
ages suggest, our results demonstrate that a degeneracy ex- 
ists between such complex structures and simple morpholo- 



gies, with VLBI observations not being able to distinguish 
between the two. 

For each element in our simulated array we also assume 
a fixed system temperature (T sya ). However, in a real ob- 
servation T sys will increase (e.g. as the antenna elevation 
decreases) and will affect the noise levels in the observations 
more than is accounted for in our simulations which are 
not prone to these effects. Note that as the source becomes 
fainter and the signal to noise level decreases, the noise in 
the image exacerbates the apparent fragmented shell. 



3.4 Residual structures 

Figure 3 shows the simulated VLBI images (presented in 
Figure 2) subtracted from the input model. To facilitate 
this comparison the input models are convolved with cir- 
cular Gaussians with sizes equal to those of the restoring 
beams used to produce the simulated images. The restoring 
beams for each epoch are listed in Table 1. The images are 
carefully aligned to give identical central pixel values and 
pixel increments. 

The residuals from the subtraction are contoured and 
superimposed over the simulated VLBI image, with contour 
levels of (-5, -4, -3, 3, 4, 5, 6, 7, 8) x a. The underlying 
greyscales are the same as those presented in Figure 2 and 
the values of a are in accordance with those presented in 
Table 1. 

By contouring the residuals above the 3<r level the struc- 
ture that is revealed can be considered to represent genuine 
artifacts arising due to the iro-plane sampling, the deconvo- 
lution process, or both. Note how the brightest features in 
the simulated image often correspond to negative regions in 
the residual map. The effects are less pronounced in the ear- 
lier epochs which exhibit higher signal to noise levels for the 
source, as well as having more compact angular sizes. Even 
in the final and faintest epoch there is residual structure at 
the 7a level. 

We quantify the level of the residual structure by mea- 
suring the rms of the signal over the region occupied by the 
shell. This is presented in Table 2 both as an absolute value 
and in units of the background rms (a). The peak and min- 
imum flux densities are also presented in order to represent 
the extremes of the image contamination. 
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Figure 2. Simulated VLBI images at nine successive epochs. Contour levels are 1%, 2%, 4%, 8%, 16%, 32%, 45%, 64% and 90% of the 
peak brightness starting with the first contour level above 3<r. The greyscale is from 8% to 100% of the peak brightness for t < 1253 
days, and from 16% to 100% for later epochs. The greyscale pixel values above the images are in units of mjy/beam, except for the final 
three epochs which are in /ijy/beam. 



4 ILLUSTRATIVE SIMULATIONS 



The first part of this section presents simplified simulated 
observations to demonstrate the coupling between the uv- 
plane sampling and the resulting image. The second subsec- 
tion takes a more realistic simulation and shows how subtle 
changes to the visibility data (e.g. antenna flagging) affect 
the image. 



4.1 The dependence of azimuthal structure on 
hour angle 

As a simplified illustration of the effects of incomplete uv- 
sampling on azimuthal structure, Figure 4 shows a simu- 
lated continuous six-hour VLBI simulation for seven differ- 
ent hour-angle ranges. Note that the tracks are truly con- 
tinuous; there are no gaps in the visibilities. These simu- 
lations are constructed with an antenna configuration and 
input model both consistent with those of day 2787 of the 
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Figure 3. Residual maps created by subtracting the simulated VLBI image (shown in Figure 2) from an appropriately convolved input 
model. Contours show the difference map, superimposed on a greyscale plot of the simulated VLBI image. Contour levels are (-5, -4, -3, 
3, 4, 5, 6, 7, 8) times the rms background level (<r) for each specific image. The greyscale runs from 8% to 100% (t < 1253 days) and 
16% to 100% (thereafter) of the peak brightness in the simulated image. This is shown in the scale above each image, the units of which 
are mjy / beam. The a values are consistent with the 3<r values presented in Table 1. 



actual VLBI observations of SN1993J. The left hand column 
shows the uu-coverage and the right hand column shows the 
resulting map. 

The most notable feature in this series of images is that 
as the uw-coverage rotates with hour angle, so does the struc- 
ture of the corresponding image. The position of the peaks on 
the ring rotate, as does the central 'bar' in the image. Note 
that the input model is purely circularly symmetric with no 



angular structure variations (see Figure 1), therefore none 
of the complex structure in the final images is real 

In order to quantify the azimuthal structure that is im- 
parted to an image, and demonstrate how it is radically 
affected by the w-coverage, the pixel values as a function of 
angle are presented in Figure 5. This simulation is carried 
out 24 times with the hour-angle coverage drifting by 30 
minutes each time, using day 2787 as the input model. The 
resulting images are 512 x 512 pixels in size, and the pixel 
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Figure 4. Simulated continuous six hour VLBI observations of the model shell corresponding to day 2787, for seven different hour-angle 
ranges. The hour-angle coverage for the top row is 00:00 - 06:00 and with each row the start and end time shift forwards by two hours. 
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Table 2. Peak and root mean square flux levels over the emitting region of the shell. 
The latter is expressed both as an absolute value and in terms of the background 
noise (cr; 3cr values are listed in Table 1.) . 



Observation 
time (days) 


Peak flux 
(mjy/bcam) 


Minimum flux 
(mjy/beam) 


Residual RMS 
(mjy/bcam) 


Residual RMS 
1" 


352 


0.392 


-0.711 


0.194 


6.7 


390 


0.937 


-0.989 


0.233 


4.9 


686 


0.073 


-0.235 


0.056 


4.4 


996 


0.316 


-0.346 


0.089 


1.3 


1107 


0.233 


-0.255 


0.073 


1.5 


1532 


0.485 


-0.527 


0.142 


1.5 


1893 


0.139 


-0.167 


0.047 


2.2 


2525 


0.108 


-0.109 


0.033 


2.7 


2787 


0.177 


-0.232 


0.061 


2.7 



values are measured around a ring with a radius of 141 pix- 
els with respect to the image centre, which on these images 
corresponds to 3.8 milliarcseconds. This value was chosen 
as it corresponds to the radius of the brightest annulus of 
emission on the input model. 

Plotting the brightness distribution in this way allows 
one to clearly see how the structure rotates with the uv 
coverage. The brightest peaks in the image appear as dark 
patches and the rotation manifests itself here as the diagonal 
striping in Figure 5. 



4.2 Subtle itti-plane changes, substantial image 
plane implications 

To illustrate how apparently subtle changes to the visibil- 
ity data can alter the structure in the final image, Figure 6 
presents uv plots and corresponding images for three differ- 
ent simulated observations. These are created using the sim- 
ulation corresponding to day 2525. Contour and greyscale 
levels, and the restoring beam, are consistent with the day 
2525 simulation presented in Figure 2. 

The left hand pair shows the simulation with the 15 
degree elevation cut off applied. For a source at the Declina- 
tion of SN1993J the effect this cut off has on the uv coverage 
compared to that of a zero degree elevation cut off is mini- 
mal; in our simulation it reduced the number of visibilities 
by less than 3%. The effect this elevation cut off has on the 
plot of uw-plane coverage is subtle and difficult to notice by 
eye. As such the complete plot is not presented here. 

To illustrate potential problems that may be caused by 
local geographical features the centre pair of images sim- 
ulation has an asymmetric (+35°/+15°) elevation cut off 
imposed on the Mauna Kea station. This drastically affects 
the longest baselines and consequently alters the structure 
in the corresponding image. For the final image the St. Croix 
station is entirely removed, which mainly affects the inter- 
mediate length baselines as shown by the two prominent 
gaps in the uv plot. 

Further simulations can be found in Appendix A where 
we present a range of simulated images corresponding to the 
day 2787 epoch for a range of Declinations, track lengths and 
hour angle ranges. 



5 CONCLUSION 

We have demonstrated in this paper how sparse uv cover- 
age can impart apparently complex azimuthal structure to 
a radio brightness distribution that is in reality morpho- 
logically simple. Rotation of the ww-plane sampling pattern 
determined by varying hour angle coverage leads to a corre- 
sponding rotation in the radio image. 

This is just one example of the kind of structures that 
are vulnerable to errors arising from sparse ww-plane sam- 
pling. For a full discussion of other structures and their vul- 
nerabilities we refer the reader to the works by Cornwell 
(1999) and Briggs (1995) who investigate the use of differ- 
ent deconvolution methods which further undermine these 
structures and cannot completely overcome the problems 
caused by incomplete information. 

By reconstructing the uw-plane distributions achieved 
in the SN1993J VLBI campaign and convolving these with 
a perfect spherically symmetric model of the supernova, we 
have shown that sparse uv coverage may be wholly or par- 
tially responsible for the evolving complex azimuthal struc- 
ture observed in real supernova remnant. 
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APPENDIX A: TRACK LENGTH AND 
DECLINATION DEPENDENCE 

The parameter space defined by factors which will affect im- 
age fidelity is enormous. Both the realistic and illustrative 
simulations presented in this paper represent a very specific 
case in terms of source Declination, track length and there- 
fore Mii-plane coverage. 

In this section we present simulations of the day 2846 
model at Declinations of 20, 40, 60 and 80 degrees with 6, 9 
and 12 hour continuous observations. Please refer to Figure 



1 to see the morphology of the input source. Seven unique 
hour-angle ranges are simulated for each declination-track 
length pair resulting in 84 unique observations. Simulated 
images are grouped by Declination and are presented in Fig- 
ures Al, A2, A3 and A4. 

The greyscale range in each image is from 16% to 100% 
of the peak value and the contours are 32%, 45%, 64% and 
90% of the peak flux. The starting value of the hour-angle 
range is listed on the left hand side of the figures, with the 
end hour angle being the start value plus the track length, 
which is printed at the top of each page. 

Imaging of the simulated data sets was automated in 
this case. Cleaning was not carried out interactively and the 
clean cycle was terminated after reaching 1000 components. 
The Gaussian fitted to the half power point of the central 
lobe of the point spread function for each image was used 
as the restoring beam in each case. This is plotted in the 
lower left of each figure. The favourable Declination (+69 
degrees) of SN1993J facilitates the use of a circular restoring 
beam as the fitted beam itself is nearly circular. For lower 
Declinations this is not applicable, as can be easily seen by 
the elongated beams of lower Declination simulations. 

Note the wide variety of image artifacts that can arise 
due to the uv-plane sampling. Very sparse coverage at low 
Declinations leads to 'polar cap'-type features, striped gaps 
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Figure 6. The upper row shows the uv plane coverage for day 2525 with various flagging conditions applied. Left to right, these are 
(a) with archival flag tables and a 15 degree antenna elevation cut off applied, (b) with asymmetric Mauna Kea elevation flagging, and 
(c) with the St. Croix antenna removed. The lower row shows the corresponding images. The contouring regime is consistent with the 
simulated VLBI images, in that the greyscale runs from 16% to 100% of the peak brightness (listed in Table 1) and the contour levels 
correspond to 1%, 2%, 4%, 8%, 16%, 32%, 45%, 64% and 90% of the peak brightness, starting with the first contour level above 3a. 
Left to right, the peak brightnesses for these images in mjy/beam are 0.65, 0.65 and 0.66. The greyscale pixel values arc in units of 
mjy/beam. 



in the uw-coverage result in dual ridge features, and central 
bars across the shell are a regular outcome. Only the 12- 
hour tracks at high Declination, featuring almost complete 
coverage, result in a ring structure with a central minimum. 
Artificial azimuthal structure is present in each of these sim- 
ulated images, even those with excellent uw-coverage. 

This paper has been produced using the Royal Astronomical 
Society /Blackwell Science KtgX style file. 
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Figure Al. Simulations of the day 2787 model at 20 degrees Declination. 
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Figure A2. Simulations of the day 2787 model at 40 degrees Declination. 
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Figure A3. Simulations of the day 2787 model at 60 degrees Declination. 
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Figure A4. Simulations of the day 2787 model at 80 degrees Declination. 



